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基于 椭 球 拟 合 的 人 体 一 服装 碰撞 检测 方法 
柳 随 , 修 进 '， 钟 李涛 ， 张 ” 娟 ， 唐 源 捍 


(西南 交通 大 学 信息 科学 与 技术 学 院 , 成 都 611756) 


摘 要 : 为 了 实现 服装 念 真 中 服装 与 人 体 的 快速 碰 挤 检测 ， 提 出 了 一 种 基于 椭 球 拟 合 的 碰撞 检测 方法 。 首 先 ， 以 测 地 
距离 等 值 线 为 基础 数据 ,结合 人 体 尺寸 与 身高 的 线性 关系 和 人 体 结 构 特征 提取 模型 特征 点 , 实现 模型 语义 分 割 ; 然后 ， 
以 径 向 距离 的 平均 值 作为 椭 球 与 模型 之 间 的 拟 合 误差 ， 采 用 前 杆 优化 的 二 分 KK 均值 聚 类 算法 逐步 增加 聚 类 中 心 个 数 
实现 对 人 体 模型 的 快速 权 闫 并 生成 一 系列 各 近 模型 的 最 小 体积 包 国 类 于 ;最 后， 使 用 生成 的 包 国 精 球 公共 人 体 楼 型 与 
布料 进行 碰撞 检测 。 实 验 结果 表明 ， 该 方法 不 仅 能 快速 实现 对 三 维 人 体 模型 的 高 度 拟 合 ， 而 且 有 效 提 高 了 碰撞 检测 的 
计算 效率 。 
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Cloth-body collision detection method based on ellipsoid fitting 


Gong Sui, Hou Jini, Zhong Litao, Zhang Juan, Tang Yuanhao 
(School of Information Science & Technology, Southwest Jiaotong University, Chengdu 611756, China) 


Abstract: To realize the rapid collision detection between cloth and human body in cloth simulation, a collision detection method 
based on ellipsoid fitting is proposed. Firstly, the geometric distance isoline was used as the basic data. And the model feature 
points was extracted by combining the linear relationship between body size and height with the structural features of the human 
body, which realized the semantic segmentation of the model. Secondly, the average of the radial distance was used to represent 
the fitting error between an ellipsoid and the model. And the number of cluster centers was gradually increased by bisecting K- 
means clustering algorithm based on pruning optimization. It could realize the rapid clustering of human model and generate a 
series of minimum ellipsoids. Finally, cloth grid points was detected with the ellipsoids instead of the model. Simulation 
experiments show that this method not only realizes the fast fitting of the human model, but also effectively improves the 
computational efficiency of collision detection. 


Key Words: collision detection; geometric distance; model segmentation; minimum volume ellipsoid; bisecting K-means 
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0。 到 诗 aligned bounding box, AABB) 层 次 结构 和 GPU 多 线程 实现 了 对 
一 高 分 辩 率 布料 的 动态 实时 仿真 。Du 等 人 句 同 样 采 用 AABB 包 


在 虚拟 服装 仿真 系统 中 ， 磁 撞 检 测 是 决定 仿真 中 视觉 效果 天 盒 处 理 多 层 服装 间 的 穿 透 ， 但 当 服装 粒度 过 大 时 依然 存在 # 
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逼真 度 的 关键 因素 ， 同 时 也 是 最 为 耗 时 的 部 分 ， 如 何 快速 检测 透 问 题 。 针 对 AABB 包围 盒 紧 密 性 差 的 问题 ，Shen 等 人 中 使 用 
出 服装 与 人 体 模型 之 间 的 碰撞 并 作出 响应 是 影响 试 衣 实 时 性 的 有 向 包围 盒 (orientedbounding box,OBB ) 和 包围 球 来 构成 双重 
一 个 关键 问题 。 为 此 ， 大 量 研究 人 员 进 行 了 深入 研究 ， 并 取得 包围 盒 层 次 结构 ， 首 先 通 过 包围 球 剔 除 不 会 相交 的 对 象 ， 再 通 
了 较 好 的 研究 成 果 上 3 。 过 紧密 性 较 好 的 OBB 包围 盒 进 行 精确 碰撞 计算 ， 该 方法 综合 


前 已 经 形成 了 两 类 比较 经 典 的 算法 : 层次 包围 合法 MI 和 ”了 OBB 的 紧密 性 和 包围 球 的 简易 性 ， 相 比 单一 的 OBB 算法 更 
空间 分 解法 。 其 中 前 者 被 广泛 应 用 于 复杂 环境 中 空间 几何 模型 ” 高 效 , 但 这 种 方法 并 不 适用 于 近 距 离 接触 物体 之 间 的 碰撞 。Guo 


之 间 的 碰撞 检测 。Salazar 等 人 外 通过 使 用 轴 向 包围 盒 (axis 等 个 则 将 AABB 拉 例 和 离散 方向 包 目 
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orientation poly-topes, K-DOP) 混 合用 于 加 速 虚拟 手术 中 器 官 形 


变 产生 的 碰撞 检测 ， 在 不 影响 精度 的 前 提 下 提高 了 碰撞 检测 效 


率 。 唐 等 人 中 利用 多 核 CPU 的 


检测 效率 ， 但 在 核心 数目 较 大 时 


行 性 ,结合 


SIMD 指令 进一步 
优化 k-DOP 模型 , 有 效 提 升 了 多 核 体系 架构 下 柔性 织物 的 碰撞 
， 会 造成 SIMD 指令 获得 的 加 


速 比 下 降 等 问题 Zhang 等 人 中 通过 使 用 多 核 CPU 和 多 核 GPU 


进行 并 行 碰撞 检测 
出 一 种 针对 分 布 式 存储 GPU 
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， 该 方法 简化 了 计算 的 了 
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法 。Ye 等 人 0 
透 现象 ， 对 于 复杂 模型 表 


交 的 模型 。 


[ 作 量 。 


Du 等 人 09 提 


的 并 行 连 续 硼 
] 采 用 基于 离散 交叉 检测 的 方法 处 理 模 型 之 间 
面 的 变形 ， 该 方法 相 比 于 传统 连续 碰 


撞 检 测 加 速算 
的 渗 


撞 检 测 方法 具有 较 好 地 性 能 ， 但 该 方法 并 不 能 有 效 处 理 
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对 于 服装 -人 体 之 间 的 冲突 检测 而 言 , 服装 与 人 体 之 间 较 为 
紧密 ， 其 运动 复杂 且 极 易 变 形 ， 采 用 层次 包围 盒 的 方式 在 运动 


过 程 中 需要 实时 更 新 层次 结构 ， 计 算 带 来 的 时 间 消 耗 成 为 上 


方法 效率 的 瓶颈 之 一 。 


近似 方法 。 王 等 人 03 提 出 一 种 三 维 


模型 球体 逼近 方法 ， 在 三 维 
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鉴于 球体 包围 盒 的 方向 无 关 性 和 碰撞 响 
应 的 简易 性 ，Wang 等 人 [3 提出 一 种 通过 球体 近似 固体 模型 
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预 处 理 阶段 
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由 于 单一 椭 球 体 难 以 实现 对 复杂 三 维 模型 的 有 


到 1 所 示 。 


原始 网 格 模型 


| 模型 分 类 | 
| 梢 球体 初 给 拟 合 
| 要 球 体 细 分 优化 
| 而 擂 检 测 与 处 理 | 
| 仿真 结果 | 
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碰 手 检测 方法 


效 拟 合 通 


党 需要 对 复杂 三 维 


模型 进行 简化 分 割 ， 使 其 


模型 体 素 化 的 基础 上 构建 内 部 


球 通 近 ， 


数量 较 多 ， 不适/ 
体 着 装 时 服装 质点 与 相连 三 角 


于 实时 性 较 高 的 碰撞 检测 。 
面 片 上 


但 该 方法 构建 的 内 部 : 


的 关系 ， 提 出 一 种 相干 碰撞 


中 
Li 等 人 0 基于 人 


检测 方法 ， 其 性 能 相 比 于 传统 


用 体 紧 密 性 差 、 


层次 结构 复杂 的 问 
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与 此 同时 ， 针 对 人 体 网 格 模型 
题 。 杜 
斯 混合 模型 聚 类 分 割 。 王 等 人 
值 边界 检测 的 三 角 网 格 模型 分 割 


盒 方 法 提升 
缺少 足够 语义 信 
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提出 一 种 基于 四 是 型 号 的 最 


算法 ， 该 方法 有 较 高 的 鲁 棒 


求 体 作 为 基 元 采用 


性 ， 但 在 不 同 场合 下 需 人 工 适 当 调 


自 适应 椭 球 包 


整 阐 值 。L 


自 项 向 下 的 方式 对 给 定 模型 构建 


于 人 体 结构 特征 ， 
寸 与 身高 有 之 间 
2.1 特征 点 提取 


模型 中 ， 椭 球体 对 模型 四 上肢 


于 等 人 PMV 提 出 的 方法 


以 测 地 距离 作为 基础 数据 ， 


+H 和 


分 解 为 一 系列 具有 
相对 较 高 拟 合 度 的 子 区 域 ， 从 而 改善 整体 拟 合 的 效果 。 在 人 体 
有 较 好 的 拟 合 性 ， 为 此 ， 本 文 基 
结合 人 体 部 位 尺 


的 线性 关系 ， 实 现 快速 有 效 的 模型 语义 分 割 。 


自动 


本 文采 | 


是 取 人 体 模型 末端 特 


征 点 ， 即 四 有 上肢 林 端 和 头顶 5 个 部 位 的 特征 点 。 此 外 ， 由 于 人 体 


文 对 GB/T 10000-88《 中 国 成 年 人 人 体 尺 寸 》 RJ 和 


92《 工 作 空 间 人 体 斥 寸 》 进 行 看 


iu 等 人 07 使 用 
层次 包围 盒 ， 


该 方法 能 有 效 逼 近 但 计算 成 本 较 大 。Meng 等 人 08 根 据 人 体 骨 
骼 结构 对 人 体 模 型 进行 初始 化 分 ， 之 后 针对 初始 化 分 
骨骼 结构 。 


的 部 位 实 
青 等 人 09 基 


羽 盒 ,然后 根据 表 


撞 响 应 ， 该 方法 


岗 椭 球 体 拟 合并 细 分 ， 但 该 方法 依赖 于 
于 改进 的 均值 聚 类 算法 构建 

面 误差 分 为 两 类 分 别 进行 如 

算 较为 耗 时 。 本 文 方法 基于 椭 球 体 对 
特性 ， 在 模型 语义 分 割 的 基础 上 进行 相 


干 和 四 上 胶 采用 不 同 的 初 


二 分 KK- 均值 聚 类 算法 实现 人 体 模型 的 


撞 检 测 ， 得 到 了 较 好 的 仿真 结果 。 
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本 文 算法 分 为 预 处 理 和 动态 模 殷 
也 距离 算法 


提取 模型 特征 点 ,采用 测 


动态 模拟 阶段 主要 分 为 模型 拟 合 币 


于 反复 El 洲 ? 村 


面 物体 
随 球 体 初 始 拟 合 ， 针 对 躯 
台 聚 类 中 心 和 大 值 ， 并 结合 剪 枝 优化 的 
R 速 椭 球 拟 合 及 布料 碰 


良好 的 表面 拟 合 


两 个 阶段 ,在 预 处 理 阶段 ， 
实现 复杂 模型 的 语义 分 割 。 


撞 响 应 两 个 方面 ， 


首先 对 


分 割 好 的 模型 数据 采用 最 小 体积 包围 椭 球 进行 初始 拟 合 ， 然 后 


针对 拟 合 误 差 较 大 的 椭 球 体 进行 细 分 优化 ， 最 后 使 / 


优化 后 的 


顶 球 体 代 蔡 原始 模型 与 布料 进行 胡 


撞 检 测 处 到 


得 到 最 终 的 仿真 


点 、 中 


AAA 


特征 
身高 
科 # 
膝 点 和 颈 点 )。 


表 1 人 体 部 分 尺寸 与 身高 之 间 的 线 | 


日 之 间 的 线性 关系 ， 如 表 1 
if 彻 (Dijkstra) 近 似 测 地 距离 算法 快速 定位 中 


究 分 析 , 将 人 体 模型 分 为 末端 
层 特征 点 和 内 层 特 征 点 三 类 , 计算 出 人 体 部 位 尺寸 与 


比例 具有 一 定 的 稳定 性 ， 不 会 随 着 身高 的 改变 而 快速 变化 ， 本 


GB/T 13547- 


所 示 。 并 以 此 为 基础 ， 结 合 迪 


屋 特 征 点 《 肘 点 、 


生 关 系 


项 目 KHJ) ”成 年 男性 (18-60 岁 ) ”成 年 女性 (19-60 岁 ) 

头 部 长 Ll 0.0836H+168.9137 0.0994H+142.6445 

前 臂 长 L7 0.5510H-283.1964 0.5672H-307.5205 

小 腿 长 L4 0.7447H-257.6398 0.7134H-207.8311 

区 高 IL3 0.8044H-243.0916 0.7683H-170.2316 

向 高 ”L2 0.9630H-311.3857 0.9309H-245.0959 

由 于 本 文采 用 的 模型 沿 y 轴 正 方向 站 立 ， 故 只 需 对 所 有 坐 
标点 按照 其 y 坐标 值 查找 出 最 大 值 > 和 最 小 值 > ，， 即 可 求 
出 人 体 身高 : 太 =|y,, 一 y,,|。 之 后 分 别 以 末端 特征 点 为 基准 点 


求解 模型 上 任意 一 点 v 与 各 末端 特征 点 p 之 间 的 测 地 距离 。 


其 中 ，y, 为 模型 
为 两 点 之 间 的 测 


G(x) = geod(v,,p,), ne(0,N) 


鉴于 人 体 模 型 是 


不 同 疏 密 程 度 的 三 4 


(1) 


的 第 nn 个 顶点 ， PP 为 某 一 末端 特征 点 ， GO 
也 距离 ， N 表示 模型 的 总 顶点 个 数 。 
网 格 构成 ， f(H) 
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所 在 截面 与 模型 相交 时 必然 与 边 相 交 ， 为 了 保证 计算 精度 ， 需 
要 求解 截面 与 三 角 网 格 的 交点 ; 检测 模型 上 所 有 三 角形 ， 求 解 
GO = 了 (五) 的 点 ，f(H) 表示 人 体 部 分 尺寸 与 身高 及 之 间 的 
线性 关系 。 如 图 2 中 AABC , 设 G(4) > f(H),G(B)<f(H); 根 
据 式 (2) 可 以 计算 出 点 d 坐标 。 同 理 , 点 q, 坐标 也 可 以 得 到 。 


-ACD 6) 人 A)+A 0) 
G(B)—G(A) 
A 
d, d, 
B C 


图 2 测 地 距离 等 值 点 


通过 上 面 的 计算 ， 可 以 获得 肘 点 、 膝 点 和 有 颈 点 所 在 的 测 
距离 等 值 线 。 会 阴 点 即 为 两 脚尖 特征 点 之 间 测 地 距离 长 度 一 
所 对 应 的 v 点 ， 模 型 的 腿 长 L5 则 用 会 阴 点 到 脚尖 点 的 测 地 昌 
离 G。( 代替 ， 即 L5 = G, (CD 。 为 了 准确 的 划分 模型 手臂， 结 
合 人 体 模型 结构 特征 ， 可 通过 点 到 直线 的 距离 准确 提取 ， 如 图 
3 所 示 ， 为 手 肘 对 应 测 地 距离 等 值 面 上 一 点 到 会 阴 点 的 最 小 
欧式 距离 点 ; v 为 以 右 脚尖 末端 特征 点 为 基准 点 ，FCz) = Z3 计 
算 的 等 值 线 上 距离 会 阴 点 欧式 距离 最 大 值 点 。 依 次 比较 六 到 
之 间 测 地 线 经 过 的 所 有 顶点 中 点 ， 到 直线 ww 的 投影 距离 4 ， 
则 投影 距离 最 大 的 v 点 即 为 腋 下 点 ，w 到 对 应 指 尖 特征 点 的 测 
地 距离 即 为 模型 手臂 长 L6= G(x)。 

2.2 ”模型 分 割 

得 到 人 体 语 义 特 征 点 以 后 ， 可 以 通过 测 地 距离 完成 三 维 人 
体 模型 的 层次 分 割 。 为 了 描述 方便 ， 设 三 维 模型 所 构成 的 顶点 
集合 为 $={v,w,…v.…svy} ， 其 中 为 模型 第 ;个 顶点 ，N 为 
三 角 片 总 数 ， 模 型 分 割 的 最 终 目 的 是 将 三 维 模型 所 有 项 点 8 按 
照 其 结构 特性 进行 划分 , 以 手臂 分 割 为 例 , 对 于 末端 特征 点 P ， 
以 及 相 邻 内 层 特 征 点 w ， 可 知 属于 手臂 部 分 的 顶点 为 
geod(v,, p;) < geod(v', p;) G3) 

类 似 地 ， 在 特征 点 的 约束 下 ， 可 以 得 到 腿 部 、 头 部 的 分 制 
结果 ， 余 下 顶点 构成 人 体 模 型 的 躯干 部 分 ， 之 后 依据 腰 高 、 胸 
高 和 左右 腋 下 点 连 线 中 垂 面 把 躯干 部 位 划分 为 6 个 部 分 。 最 后 
针对 已 经 分 割 好 的 模型 四 肢 ， 找 到 相应 的 中 层 特征 点 ， 形 成 一 
个 需要 二 次 分 割 的 子 部 分 集合 , 根据 式 (3 ) 再 次 分 割 为 两 部 分 ， 
完成 模型 的 语义 分 割 ， 由 图 4 可 以 看 出 ， 男 士 模型 和 女士 模型 
的 分 割 结 果 都 能 能 较 好 的 反映 人 体 结构 特征 。 


3 ” 椭 球 体 拟 合 及 优化 
3.1 最 小 体积 包围 椭 球 
鉴于 AABB 包围 盒 对 曲面 模型 拟 合 度 较 低 和 OBB 包围 盒 


碰撞 检测 相交 测试 较 多 的 特性 ， 本 文采 用 对 弧度 表面 更 具 拟 合 
性 的 椭 球 体 作为 模型 的 基 元 包围 体 。 采 用 文献 [19] 中 的 方法 对 


Es 


IE 长 


上 


Sr 


三 维 模型 点 集 s 生成 一 系列 最 小 体积 包围 


iv 合作 期 刊 
装 A 检测 


椭 球 (minimum 


volume enclosing ellipsoids, MVEE), 
即 MVEE(S)， 如 下 性 质 : 


在 其 


基础 上 进行 改进 ， 


MVEE(S)=s0. ={xeR’:(x-c) Q(x-c)<1 (4) 


MVEE(S)/d Sconv(S)c MVEE(S) (5) 


式 (4) 中 ， so, 表 示 椭 球 ， 由 3x3 维 


来 表示 ， 其 中 的 特征 向 量 和 特征 值 分 别 为 


和 半 长 轴 长 度 。 式 (5) 中 : MYVYEE(S)/d 表示 围 


的 正定 矩阵 g 和 中 心 c 


随 球 的 半 长 轴 方 向 


绕 MVEE(S) 的 中 


心 对 其 缩放 1/d ，comyv(8) 是 $ 的 凸 包 。 


经 过 预 处 理 的 人 体 模型 被 按照 语义 分 为 多 个 不 同 的 区 域 
PS={5,,S,,…s5,,…S,}， 保 证 了 分 割 结果 和 人 体 拓扑 结构 的 一 致 


z 


= 


图 3 人 体 尺 寸 及 会 阴 点 、 腋 下 点 示意 图 


个 个 


加 4 ”模型 语义 分 割 及 会 阴 点 、 


3.2” 椭 球体 拟 合 误 差 
在 对 子 区 域 $ ,的 椭 球 体 拟 合 过 程 


行 数值 化 表示 。Simari 等 人 3] 通 过 使 月 


对 应 的 权重 ,实现 了 较 好 的 拟 合 效 果 。 
与 椭 球 体 体积 之 间 的 误差 作为 自 适应 成 


| ， 为 了 有 效 评估 生成 本 
球体 的 拟 合 情 况 ， 需 要 对 拟 合 椭 球 与 子 区 域 之 间 的 近似 误差 进 


披 下 点 提取 结果 


有 欧式 径 向 距离 、 表 面 法 
线 和 曲率 的 加 权 平 均值 作为 拟 合 情 况 的 判断 标准 
Lu 等 人 P 基 于 原始 模型 


通过 控制 相 


分 三 


维 模型 ， 生 成 紧密 


包围 的 椭 球 , 但 计算 模型 的 体积 较为 复杂 , 计算 消耗 太 大 。 由 于 
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本 文中 已 经 对 人 体 模型 语义 分 割 处 理 ， 椭 球体 已 具备 
局 部 区 域 椭 球 进行 
值 。 如 式 (6) 所 示 ， 拟 合 误差 可 通 


合 特性 , 仅 对 


与 近似 表 


之 间 径 向 距离 的 平均 


其 中 : E, 表示 模型 子 
面 点 中 第 :个 顶点 ， 开 为 子 
田 上 和 


示 顶 点 ， 在 李 


E, 


球体 表 


1 K 
2 


3.3“” 椭 球体 细 分 优化 


人 体 模型 
合 表面 特征 的 


经 过 语义 分 割 和 初始 拟 合 之 


田 化 处 到 
过 计算 人 体 模 型 子 区 
值 来 近似 : 


的 径 向 投影 点 ， 


Vi 9s 0 


椭 球 ， 但 是 在 模型 役 干 部 


yy 


大 的 椭 球 体 ， 


模型 的 过 度 下 


雹 化 把 


6 个 子 区 域 中 
剪 枝 优 
的 椭 球 


定 值 为 止 。 


这 主要 是 


的 精确 性 ， 大 大 降 1 


氏 了 碰撞 检测 的 效率 。 


良好 的 拟 
E, 无须 准 确 地 求解 误差 
或 表面 点 


(6) 


区 域 8 的 椭 球 体 拟 合 误差 v, 为 子 区 域 表 
区 域 中 包含 的 顶点 总 数 ， 9 (v) 表 
| 表示 欧 氏 距离 。 


后 ， 得 到 了 一 系列 符 
依然 存在 拟 合 误差 较 
于 在 模型 服 干 部 位 细 分 级 数 不 够 所 致 ， 
可 以 通过 增加 全 局 和 迭代 次 数 来 减 小 拟 合 误差 ， 但 同时 
分 而 引入 大 量 的 拟 合 椭 球 


也 会 造成 


体 ， 同 时 不 能 保证 拟 合 


为 了 实现 躯干 全 


可 琶 


合 ， 本 文 依据 腰 高 、 胸 高 和 左右 
干部 位 划分 为 6 个 部 分 ，K 均值 聚 类 的 初始 磊 值 设 为 6， 相 应 


腋 下 点 连 线 中 于 


E 面 把 抠 


包含 的 顶点 的 均值 坐标 作为 初始 质心 ， 然 后 基于 


此 的 均值 聚 类 算法 


体 结构 。 剪 枝 优化 算法 思想 如 下 : 
态 下 ， 所 有 的 点 都 不 会 从 一 个 角 转 移 至 
生 改 变 ， 即 : K 均值 最 明显 的 收敛 过 程 会 发 和 9 


进行 二 次 划分 ， 


得 到 逼近 效果 较 好 
由 于 K 均值 在 收敛 状 
| 另 一 个 戏 ， 质 心 不 再 发 
在 算法 运行 的 前 


期 阶段 ， 故 在 某 些 情况 下 为 了 增加 算法 的 执行 效率 ， 可 以 改变 
原始 算法 的 结束 条 件 ， 采 用 


量 小 于 特 


个 聚 类 中 心 的 均值 变化 


如 图 


5 所 示 ， 


vs 和 wv 分别 为 子 区 域 


Yp,3 


其 中 : 丈 、 如 分 别 为 Avm 


为 预 设 阔 值 。 


为 了 保证 二 次 划分 后 的 结果 和 效率 ， 本 文采 / 


聚 类 终止 条 件 示 意图 


恨 设 区 均值 聚 类 的 初始 k 
中 各 特征 点 域 的 均值 华 标 ,vw 、v ,和 
' , 分别 为 下 一 步 对 应 的 聚 类 中 心 ， 则 聚 类 的 限定 条 件 为 


值 等 于 3 


和 Vp 、 


[5 -7l|so (7) 
52V53 和 Ay ivyv3 的 重心 坐标 ，o 
欧式 距离 作 

判断 


为 聚 类 准则 ， 在 剪 枝 优化 的 同时 逐个 增加 值 和 中 心 ， 


椭 球 体 的 拟 合 误差 玉 《是 否 满 足 ，/ 为 预 设 阔 值 ， 满 足 条 件 


则 停止 细 分 。 
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算法 的 具体 步骤 如 下 : 

a) 判 断 是 否 为 肛 干 ， 若 为 服 干 ， 则 初始 化 值 为 6， 否则 
等 于 1， 进 行 剪 枝 优化 K 均值 聚 类 ， 对 聚 类 形成 的 所 有 秘 包 含 
顶点 计算 最 小 体积 包围 椭 球 。 


b) 找 出 拟 合 误 差 最 大 的 包围 椭 球 ， 使 用 最 长 轴 的 两 个 端点 


作为 新 的 中 心 点 将 包围 椭 球 对 应 的 点 集 进行 划分 
值 增加 1。 


点 集 的 原 有 中 心 点 ,， 


， 同 时 握 弃 该 


0o) 根 据 新 的 值 和 新 的 中 心 点 对 子 区 域 中 所 有 点 进行 剪 梳 


优化 -均值 聚 类 。 
d) 对 聚 类 生成 的 k 个 点 集 重新 计算 得 到 新 的 最 小 体积 包围 


顶 球 。 


6) 比较 所 有 椭 球 包 硬 


小 ， 着 E, > 


上 述 过 程 


盒 的 殷 
人 A 则 转 到 b)， 反 之 停止 细 分 。 


点 进行 重新 归 类 时 需 添加 限定 条 件 : 


其 中 : 


Bs 


点 当前 所 在 椭 球 ， 椭 球 j 为 


4 


服装 布料 作为 一 种 高 度 可 变形 
开 紧 密 宪 盖 在 人 体 模型 表面 ， 
地 检测 才能 得 到 真实 的 效果 。 男 一 方面 ， 
同 ， 在 人 体 模 型 的 不 同 
此 ， 为 了 提高 碰撞 检测 的 效率 ， 保 证 计算 速度 和 视觉 效果 ， 


性 ， 


L; > 天 


合 误差 , 与 预 设 域 值 4 的 大 


Ph, 为 了 使 拟 合 误差 E， 收敛 ,在 聚 类 过 程 中 对 


(8) 


E, 分 别 为 椭 球 i、j(i zz j) 的 拟 合 误差 。 椭 球 ; 为 该 


其 他 的 顶 球 包围 


合 
Lo 


碰撞 检测 与 处 理 


本， 具 
因 出 


轻柔 灵活 的 材质 特 
与 模型 的 碰撞 


需要 精 ; 


二 


服装 类 型 的 不 


考虑 到 
部 位 所 需 的 精准 程度 定然 


也 不 一 样 。 因 
本 


文采 / 
拟 合 误差 ， 


划分 为 一 系列 高 度 拟 合 的 椭 球 体 ， 使 用 这 些 椭 


j 双 重 拟 合 误差 标准 ， 针 对 模型 的 躯干 和 


0 设置 不 


Hl 


的 阐 值 4 。 


于 


Ea 


型 与 服装 进行 碰 措 


能 够 有 效 提高 碰 


重 ， 


保证 仿真 的 逼真 度 。 


每 陋 一 段 时 间 间 隔 检测 碰撞 是 否 发 生 ， 一 旦 检测 到 碰撞 ， 


0 肢 采用 
基于 以 上 方法 


撞 检 测 的 效率 ， 


不 同 的 
， 人 体 模型 被 
本 取代 原始 模 
同时 也 能 


则 改变 相应 质点 的 速度 和 位 置 来 避免 穿 透 。 对 于 服装 模型 上 的 


一 点 vy 带 入 式 (9) 计 算 < ， 式 中 0 为 顶 球 的 正定 矩 
C >1 表 示 没 有 碰撞 发 生 ，“ <1 表示 发 生 碰撞 ， 点 p 


中 心 ， 当 


处 于 椭 球 内 部 ， 需 要 将 其 移动 至 距离 ] 


阵 ，c 为 椭 球 


页 点 当前 位 


最 近 距 离 的 


椭 球 体 表 面 ， 


和 顶点 当前 位 置 连 线 移动 至 椭 球 体 表 面 。 


通常 比较 简单 的 处 理 


方法 是 直接 沿 着 椭 球 体 中 心 


“=(z-c) Q(x-o) 


但 上 酝 


[的 方法 只 针对 


一 个 椭 球 体 包围 


盒 成 立 


， 当 处 理 


(9) 


多 个 


相交 椭 球 体 的 碰撞 检测 时 ， 相 应 顶点 可 能 被 移动 到 错误 位 置 。 


如 图 6a 所 示 ， 


当 服装 模型 上 一 点 p 与 F 


测 响应 时 发 4 


情况 1 顶点 从 p, 到 p 时 ，p, 处 于 多 个 机 球 相交 处 ， 


每 个 椭 球 都 检测 到 碰撞 响应 ， 不 同 的 响应 顺序 将 导致 不 同 的 响 


三 种 错误 情况 : 


相交 椭 球 进行 碰 


撞 检 


于 


应 结果 ， 易 造成 系统 不 稳定 。 
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情况 2 顶点 从 p, 到 p, 时 ，p, 被 移动 到 椭 球 0, 的 表面 ， 
于 是 该 点 又 被 移动 到 椭 球 0, 的 


但 是 仍然 处 于 椭 球 0, 的 内 部 ， 


表面 ， 然 而 因 
球 0 内 部 。 


情况 3 顶点 从 p, 到 p, 时 ,顶点 被 移动 到 椭 球 错误 的 一 侧 。 


Va Ey 
h <.> 
/~ 


So 


(a) 错 误 位 置 


为 O, 部 分 存在 了 


F oO 内部， 最 后 该 点 依然 处 了 


nm 


李 


图 6 


(b) 位 置 修正 


针对 上 面 可 能 出 现 的 问 


可 起 ， 


首先 找 出 所 有 与 该 点 发 生 和 而 


撞 


点 位 


a 


前 点 与 该 点 初 


6 (b) 所 示 ,， 和 若 质点 
上 磁 撞 ， 
修正 后 的 位 置 ， 


六 =V—2:(v.n,): 


其 中 : 
球 心 0, 到 Pp 的 向 量 。 


5 ”实验 结果 及 分 析 


;然后 计算 这 些 相交 点 
6 位 置 距离 最 小 
球 ， 所 以 可 以 把 与 距离 最 小 椭 
pp 移动 于 
分 别 相交 于 Pp 和 p ， 根 : 


其 速度 计算 如 式 (9) 所 示 。 


的 椭 球 ， 


v. 表示 修正 后 的 速度 ，v 是 修 ] 


本 文 将 碰撞 处 理 


与 该 点 初始 位 
的 椭 球 意味 着 最 早 发 生 碰 撞 的 李 
球 的 碰撞 相交 点 作为 响应 位 置 
上 p, 过 程 中 与 椭 球 O 和 


童 检测 处 理 


分 为 两 个 阶段 ， 
计算 出 每 一 个 碰撞 相交 
之 间 的 距离 ， 磁 
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法 由 于 采用 反复 迭代 的 方 
异 较 大 ， 效 率 相对 较 低 。 


Si 


式 ， 对 于 不 


同 姿势 的 模型 拟 合 耗 时 差 


表 2 
可 以 看 出 对 


于 同一 模型 ， 


球体 代替 角色 模型 完成 与 
力 和 风力 的 作用 
色 模 型 之 间 的 相互 穿 
表 2 


耗 时 较 少 , 椭 球 拟 合 更 接近 人 体 模 型 。 


透 , 相 


所 示 为 本 文 方法 与 文献 [19] 两 种 方法 的 数据 对 比 。 
拟 合 椭 球 数 相 近 的 情况 下 ， 本 文 方法 


图 9 和 图 


9 和 图 10 展示 


服装 网 格 2 


间 的 碰撞 处 理 ， 服 装 在 重 


紧密 依附 于 人 体 模型 ， 


比 于 传统 包 硬 


eva 


d 
川 


本 文 方法 和 文献 [19] 两 种 方法 对 比 


模型 ” 顶点 数 


椭 球 数 


耗 时 /s 


文献 [19] 


本 文 (分 割 + 拟 合 ) 


4065 
3160 
4210 
5000 
8055 
7133 


5000 


60 
56 
90 
50 
59 
50 
70 


7.6 
下 


6.7 


14.2 


11.8 
10.1 
10.0 


0.8 + 1.4 
0 和 手 


0.8 十 


本 5 
寺 2.9 


十 2.6 


本 文 算法 中 将 人 体 模 


O 发 


(9) 


站 


E 前 的 原始 速度 ，n, 表示 


本 文 实 验 是 在 Window7，Intel(R) Core(TM) i3-2350M CPU 


@2.3GHz,4 GB 内 存 笔记 本 


电脑 环境 下 , 以 C++ 语言 编写 代码 ， 


不 同 的 初始 K 值 和 拟 合 
的 顶点 均值 坐标 作为 初始 
的 均值 聚 类 算法 进 和 


sx 


K 均值 聚 类 过 程 的 迭代 次 


I 


型 分 为 驱 必 


和 四 


两 大 部 分 ， 分 别 使 


合 误 差 阐 值 ， 
质心 ， 对 所 


了 模型 椭 球 体 的 
善 了 人 体 模型 躯干 部 位 椭 球 拟 合 误差 较 大 的 问题 ， 


把 躯干 各 子 
割 区 域 采 
优化 细 分 。 该 方法 有 


区 域 中 包含 
j 剪 枝 优化 
效 改 
同时 减少 了 


数 ， 有 效 提 


高 了 仿 


速度 。 实 验 中 针 


对 不 同 姿势 的 模型 , 设置 相同 的 拟 合 误差 闪 


通过 控制 参数 的 大 小 ， 分 
拟 合 效果 ， 其 中 a、 b、< 


型 ， 


相同 拟 合 阔 值 效果 ， 胃 


以 明显 的 看 到 , 在 没有 考虑 拟 合 


躯干 的 拟 合 效果 较 差 ， 细 


别 得 到 图 8 


为 男士 2 模型 ， 
左 至 右 分 别 对 应 单一 椭 球 初始 拟 合 效 果 ， 物 干 和 四 歧 设 
友 设 置 不 同 拟 合 
净值 的 初始 拟 合 

后 的 椭 球 拟 合 效 果 较 好 ， 


干 和 四 


分 优化 


值 和 聚 类 截止 闹 值 ， 
中 a、b、c、d、e、f 所 示 
d、e、 了 为 女士 1 模 


闵 值 效果 。 可 
效果 中 ， 模 型 
但 是 


编译 环境 为 Visual Studio 2010， 运 用 OpenGL 图 形 库 对 所 有 模 于 使 用 同一 拟 合 误 差 阐 值 ， 导 致 椭 球 数 较 多 ， 针 对 不 同 部 位 
型 进行 泻 染 。 本 文 使 用 的 人 体 模型 均 从 Marvelous Designer3 服 。 设置 不 同 拟 合 误差 闹 值 后 不 但 能 保证 拟 合 的 精度 ， 同 时 降低 了 
装 仿真 软件 导出 ， 服 装 布料 的 动力 学 模型 采用 文献 [25] 中 方法 椭 球 体 的 数量 ， 有 效 改善 了 碰撞 检测 的 效率 。 
构建 ， 然 后 在 重力 和 风力 P9 的 作用 下 进行 动态 仿真 。 为 了 达到 精确 碰撞 检测 的 要 求 ， 通 过 针对 不 同 的 部 位 设置 
本 文 方法 采用 人 体 部 分 尺寸 与 身高 开 的 线性 关系 提取 中 层 。 恰当 的 拟 合 误 差 病 值 ， 可 以 有 效 改善 不 同 区 域 的 精度 ， 实 现 多 
特征 点 ， 结 合 人 体 结构 特征 提取 腋 下 点 和 会 阴 点 ， 最 后 使 用 测 ，” 阔 值 拟 合 效果 ， 满 足 不同 服 装 碰撞 检测 精度 的 需求 ， 例 如 图 9 
地 距离 进行 模型 分 割 ， 尽 管 逼近 人 体 扩 寸 有 计算 误差 ， 但 是 基 ”中 宽松 裙 装 仿真 时 ， 腿 部 区 域 可 以 设置 较 大 的 拟 合 误差 阔 值 ， 
本 上 符合 人 体 拓扑 结构 ， 为 进一步 验证 分 割 方法 的 有 效 性 和 拟 。 ”相反 图 10 中 紧身 服装 试 衣 时 可 设置 较 小 拟 合 误差 闵 值 。 如 图 8 
合 效果 的 优越 性 ， 将 本 文 算法 与 文献 [19] 的 拟 合 结果 进行 实验 ”所 示 ， 通 过 对 各 部 位 设置 不 同 阔 值 标准 ， 模 型 不 同 部 位 的 拟 合 
比较 ， 采 用 不 同性 别 、 体 型 和 姿势 的 模型 分 别 进行 实验 ， 如 图 ”效果 明显 改善 ， 得 到 合理 的 划分 和 高 度 拟 合 的 椭 球 。 
7 所 示 ， 可 以 看 出 ， 两 种 方法 均 能 实现 对 不 同 姿势 的 人 体 模型 在 完成 了 人 体 模型 的 椭 球 体 拟 合 之 后 ， 本 文 使 用 高 度 拟 合 
的 高 度 拟 合 ， 但 是 本 文 方法 无 论 是 对 正常 体型 还 是 肥胖 体型 ， 的 椭 球 体 代 蔡 原始 模型 与 服装 在 风力 作用 下 的 动态 碰撞 检测 进 
分 割 结果 均 能 较 好 地 反映 人 体 层 次 结构 特征 ， 均 能 以 较 快 的 速 ” 行 了 实验 。 分 别 对 两 种 姿势 的 人 体 模型 着 装 衣 裙 进行 仿真 ， 仿 


度 完成 拟 口 ， 


且 拟 合 效果 更 佳 , 提高 了 仿真 效率 ; 而 文献 [19] 方 


真 结果 如 图 9 所 示 ， 在 重力 和 风力 作 


] 下 ， 均 无 穿 透 现象 ， 念 
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下 


真 效果 逼真 。 


a. 文 献 [19] 方 法 : 


b. 本 文 方法 : 


(d) 初 始 拟 合 (9) 同 
图 8 不 同 模 型 初始 


为 了 


10 所 示 为 身 着 紧身 衣服 的 女 了 


一 闵 值 4 


验证 本 文 方法 在 碰撞 检测 处 理 上 的 优越 性 


上 模型 和 


士 模 


型 的 仿 
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区 随 ， 等 : 基于 椭 球 拟 合 的 人 体 一 服装 碰 挤 检 


吉 疆 时 


(2) 男士 1 


(DD 多 闵 值 4 
以 合 及 细 分 优化 效果 对 比 


法 和 文献 [27] 中 算法 进行 实验 比较 , 为 了 比较 的 公司 


区 


， 将 本 文 算 


性 , 针对 相 


同 的 实验 模型 女士 2， 分 别 采 用 增加 人 体 模型 顶点 数 和 增加 服 


将 不 同 顶 点 数 及 三 角 玫 


装 模 型 顶点 数 的 方式 逐步 增加 整个 实验 着 装 模 型 的 总 顶点 数 ， 


i 片 数 的 人 体 着 装 模 型 应 用 于 文献 [27] 


方法 和 本 文 方法 ， 通 过 比较 


均 每 帧 的 


撞 检 测 处 


时 间 得 到 


图 11 和 
仅 增 加 人 体 模 型 顶点 
文 方法 每 
的 方法 呈 
段 本 文 方法 仅仅 只 需 判 


- 口 仆 。 


(3) 肥胖 女 


7 不 同方 法 覃 球体 拟 合 效果 


图 12 所 示 曲 线 对 比 图 ， 


项 碰撞 检测 处 理 
现 线性 递增 


的 时 


日 


更 新 服装 顶点 位 置 ， 省 
本 文 方法 对 相同 模型 
拟 合 椭 球 数 变化 较 小 ， 从 而 每 帧 春 
数 增加 而 增加 而 文献 [27] 将 
到 了 碰撞 检测 计算 的 过 程 中 ， 


的 顶 


存 空 间 , 但 是 也 增加 了 
出 在 仅 增 
两 种 方法 在 每 帧 确 
耗 时 较 文献 [27] 降 低 
测 的 平均 / 


6 


人 体 部 分 斥 寸 与 身高 联 的 线 局 


| 


采 


了 拟 合 效率 及 碰撞 检测 处 理 效率 ， 但 针对 不 同 模型 ， 
的 闵 值 大 小 4 需 人 为 设 定 。 除 此 之 外 ， 本 文 方法 
静态 模型 的 服装 碰撞 检测 ， 对 了 
椭 球 与 皮肤 同时 绑 定 骨骼 , 仅 对 关节 区 域 动态 更 新 椭 球 包 


月 


0 服装 模型 顶点 数 时 ， 


的 趋势 逐渐 增加 ， 原 因 在 于 而 
断 服装 顶点 与 拟 合 椭 球 的 位 
去 了 动态 更 新 包围 盒 的 时 
点 数 不 敏感 , 随 着 模型 顶 


邻接 三 角 


碰撞 检测 


AN 
[mm 


(4) 肥胖 男 


单位 为 秒 。 从 图 11 可 以 看 出 
寸 ， 随 着 着 装 模 型 顶点 总 数 的 增加 ， 本 
十 间 基 本 保持 不 变 ， 相 反 文献 [27] 


在 


撞 检 测 处 理 时 间 
对 的 基 


撞 检 测 处 理 
关系 ， 
司 。 此 外 ， 
点 数 的 增加 ， 
不 随 顶 点 
本 图 元 标记 放 


虽然 减少 了 存储 该 标记 所 用 
| 算 的 复杂 度 。 从 图 12 可 


随 着 着 装 模 型 顶点 总 数 的 


的 内 
以 看 
兽 加 ， 


时 要 


结束 语 


本 文 针对 服装 仿真 中 人 体 模型 包 
效率 问题 ， 提 出 一 种 基 了 
传统 包 


四 


F 李 球体 的 人 


撞 检测 处 理 的 时 间 都 稳定 增长 ， 但 本 文 方法 


10 倍 。 综 合 之 ， 本 文 方法 在 处 理 碰撞 检 


多 于 文献 [27] 方 法 。 


围 盒 的 紧密 性 和 碰撞 检测 


体 服装 碰撞 检测 方法 ， 与 
盒 相 比 ， 本 文 方法 在 不 需要 精确 分 割 的 前 提 下 ， 


通过 


前 枝 优 化 二 分 K 均值 欠 代 的 方式 生成 系列 紧密 椭 球 ， 


关系 ,有 效 改善 了 模型 分 割 速率 ， 


提高 


拟 合 误差 


目前 只 适 / 


于 
| 可 将 


合 


国 


» 


| 
| 

| 
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(0) 重 力 和 风力 


(Oo) 正面 (重力 ) (dd) 侧面 (重力 ) 
图 9 宽松 裙 装 在 重力 和 风力 作用 下 效果 ( 左 : 女 士 2, 右 :女士 3) 


(gd) 重 力 和 风力 
10 紧身 服装 试 衣 效果 
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